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Based on a phenomenological model with s± or s-wave pairing symmetry, the mixed-state effect on 
the low-energy spin dynamics in optimally-doped iron pnictide superconductors is studied by solving 
Bogoliubov-de Gennes equations. Our results of the spin susceptibility at q = Q in the normal, 
superconducting and mixed states agree qualitatively with recent neutron scattering experiments. 
We also propose that the field-induced intensity change shows different behaviors between the s± 
and s-wave symmetries in both momentum and real space, thus it can be used to distinguish these 
two pairing symmetries. 
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Introduction. — The superconducting (SC) pairing 
mechanism and symmetry in iron pnictide superconduc- 
tors have long been debated even since their discovery 
in 2008^ Theoretically, one commonly believed idea is 
that the Cooper pairs are formed via spin fluctuations, 
leading to the s± pairing symmetry where the paring or- 
der parameter (OP) changes sign between the electron- 
and hole-pockets^— Later, orbital-fluctuation-mcdiatcd 
pairing was also proposed and the symmetry in this 
case is s-wave without sign reversal*^ Experimentally, 
the evidences for the paring symmetry (s± or s) are 
not definite cither. For example, in optimally-doped 
Bao.6Ko.4Fe2As2 and BaFei.85Coo.i5As2, the SC gaps 
measured by angle-resolved photocmission spectroscopy 
(ARPES) are almost isotropic on both the electron- and 
hole-pocketSi 10 ^ 1 However, due to the insensitivity of 
ARPES to the sign of the OP, the paring symmetry can 
be either s± or s-wave. Neutron scattering (NS) experi- 
ments observed a resonance at Q = (±7r, 0) and (0, ±7r) 
(in the lFe/cell Brillouin zone) in the SC state^ - — 
which was initially interpreted as the evidence for the s± 
symmetry* 1 ^ - — Later it was found that the resonance can 
also emerge if the symmetry is s-wavei 22 ' 23 Furthermore, 
the local density of states probed by scanning tunneling 
microscopy (STM) with or without the magnetic field 2 ^ 
can be reproduced by assuming cither the s± or s-wave 
symmetry^ 

Recently Hanaguri et al. used STM to image the quasi- 
particle interference (QPI) and stated that the QPI in 
the presence of a magnetic field can be used to distin- 
guish the s± and s-wave symmetries^ which has been 
verified by our previous worki 2 ^ In addition, some NS 
experiments studied the effect of a magnetic field on the 
low-energy spin dynamics^ - — Specifically in Ref. H^, a 
perpendicular magnetic field applied in optimally-doped 
BaFe1.9Nio.1As2 reduces the intensity and energy of the 
resonance while broadens its width. A natural question 
arises: Is the evolution of the resonance with the mag- 



netic field consistent with the s± or s-wave scenario and 
arc there differences between them as in Ref. l26l ? In this 
paper we address this issue and propose a method to test 
it. The mix-state effect on the spin dynamics is studied 
by solving the Bogoliubov-de Gennes (BdG) equations in 
the presence of a perpendicular magnetic field and then 
conducting the real-space random-phase approximation 
(RPA) to calculate the spin susceptibility \- Previously, 
for the iron pnictides without the magnetic field, x nas 
been calculated in momentum space.— ~— However in the 
presence of a magnetic field, the translational symmetry 
is broken and \ nas to be calculated in real space which 
requires much more computational time and resources. 
Therefore a theoretical investigation of the mixed-state 
effect on the spin dynamics in the iron pnictides is still 
lacking and needs to be done urgently. By studying this, 
we find that while the spin excitation spectra in the s± 
and s-wave cases are similar to each other no matter the 
magnetic field is present or not, there are clear differences 
in the field-induced intensity change which can be used 
to distinguish the two pairing symmetries. 

Method. — We adopt an effective two-orbital model on 
a two-dimensional N xN lattice which captures the basic 
Fermi surface structures of the iron pnictides— with a 
phenomenological form for the intraorbital pairing terms. 
The Hamiltonian can be written as~r— 

H = ~ t 'ij,al3 c io a , c jP-> + '52[- V + U ( n j0p) 

+ (U - 2J H )(n j V ) +(U- 'AJH)(n jBv )]n ]Pu 
+ Y,(A ijia pcl a ^ m + H.c). (1) 

ij,a/3 

Here i,j are the site indices, a, (3 = 1,2 are the orbital 
indices, v represents the spin, fx is the chemical potential 
and rijp v = Cjp^Cjpi, is the number operator. U and 
Jh ar e the onsite intraorbital Hubbard repulsion and 
Hund's coupling, respectively. Here we have the interor- 
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bital Coulomb interaction U' — U — 2Jh according to 
symmetry.^ Ay )Qi g = Vi l*"P ((cj^Ci a f) - {cjf}^c iai )) is 
the intraorbital spin-singlet bond OP, where Vij is the 
next-nearest-neighbor [i = j ± (x ± y)] or onsite [i = j] 
attraction we choose to achieve the s± or s-wave sym- 
metry, respectively. In the presence of a magnetic field 
B perpendicular to the plane, the hopping integral is 
t 



tij,a/3£xp[i-§^ Jj A(r) • dr], where $ = hc/2e is 
the SC flux quantum, and A = {—By, 0, 0) is the vector 
potential in the Landau gauge. Following Ref. UmL 
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'-h + '-^h 



at = P,i = j±x(y), 
n = f3,i = j ± (x + y), 

\ ^ v 2 ^ ^ + ± - 1 2 ±J -h a = 0,i = j ± (x - y), 
t 4 a f3,i = j ± (x ±y), 

otherwise. 

(2) 

Equation ([I} can be diagonalizcd by solv- 
ing the BdG equations H = C'MC, where 
Cft = {■■■ , ct lf , Cjij., cj 2t , Cj-24,, • ■ • ), subject to the self- 
consistency conditions: (n^) = £) fc=1 |Qm-ifc| 2 /(-Efc), 
i/94.) = 1 — J2k \Qmk\ 2 f(Ek) and \j, a /3 = 



[nm) = 1 - Efe \Q m k?f{E k ) 

¥f S c*pY,k(QmkQnk + Q* n+ ikQm-ik)f{E k ). Here L = 
AN 2 , m = 4(j y + Nj x ) + 2f3, n = 4(i y + Ni x ) + 2a-l and 
Q is a unitary matrix that satisfies (Q^MQ)ki = SkiEk- 
Here we used i = (i x ,i y ) and j = (j x ,j y ), with 
ix, jxiiyi jy — 0, 1, . . . , N — 1. /i is determined by the 
doping concentration x through J^jpv n jpv = 2 + x. 



The bare spin susceptibility i 
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~(( C i8i C io't\ C jit C 30i) )u+ir) 
L 

= QpkQnl{Q*olQmk ~ QokQml) 

fc,i=l 

/(£*)+ /(£Q-1 

uj-Ek-Ei+iri' { ' 

where o = 4(j y + iVj'a.) + 27 - 1 and p = 4(i v + Ni x ) + 26. 
Including interactions within RPA, the full susceptibility 
can be calculated through 



x?-(t,i,w) + 5Z5^x? ? (»,«,w) 

Uh{v,v)x»Av,j,u). (4) 



Here we omitted the superscript — K The interaction 
vertex C// t is nonzero only when i = j and satisfies 

\]l = U, Ul = J H , Uh=U- 2J H , Ul = J H , (5) 

where r ^ t. The experimentally measured susceptibility 
is proportional to x (q; w )j where 



■J) 



- j,w)cosq- (i-j), 



(6) 



with f(i-j,w) = -fa ]TV Q7 (« - j, j, w) being the 

average real-space spin-spin correlation function. 

The parameters are chosen as <i_4 = 1,0.4,-2,0.04. 
N = 32 and magnetic unit cells are introduced where 
each unit cell accommodates two SC flux quanta. Vij 
[i = j ± (x ± y)] and Vu are chosen to be —2 and 
—3.2, respectively. We fix x = 0.1, corresponding to 
the optimally-doped BaFc1.9Nio.1As2, which enables us 
to compare the results directly with those in Ref. |33|. The 
temperature T is set to be and 0.7Ao- Here Ao is the 
SC OP we solved at (B = 0,T = 0). At T = 0.7A , 
the SC OP is 0, indicating that the system is in the 
normal state. In the following we calculate \ a t 
(B = 0, T = 0.7A ), (B = 0, T = 0) and (B ^ 0, T = 0), 
corresponding to the normal, SC and mixed states, re- 
spectively. 

Results. — For the s± pairing, we choose U and Jh to 
be 2 and 0.2C7, respectively. With this choice of param- 
eters, the system stays paramagnetic no matter it is in 
the normal, SC or mixed state. Meanwhile x (Q 5 W ) m 
the SC state is enhanced around uj 2Ao as compared 
to its normal-state counterpart. Both the disappearance 
of the magnetic order and the enhancement of spin sus- 
ceptibility around 2Ao are typical characteristics of the 
optimally-doped iron pnictides^ -17 ' 40 ' 41 and the results 
shown below remain valid as long as U and Jh are in 
a reasonable regime. In this case x (q> w ) m the SC 
state strongly peaks at q = Q and uj — 1.8A (uj res ), 
so in the following we normalize x (q> w ) by the value of 
X (Q, tJres) in the SC state. The constant-q scan in Fig. 
[Ha) shows that, the normal-state x (Qi w ) ( re d dash 
dot) exhibits clear gapless continuum of excitation. Upon 
entering the SC state, a spin gap gradually opens below 
to « Ao where the value of x (Q^) becomes smaller 
than that of the normal state and the low-energy spec- 
tral weight is transferred into a resonance at uj = uj res 
(black solid), consistent with previous momentum-space 
calculations without the magnetic fieldf^— thus justi- 
fying the validity of our method. While in the normal- 
state, x (q> w ) is almost unaffected by the magnetic field, 
in the mixed state, vortices emerge and the intensity 
of the resonance in x (Qj w ) is clearly suppressed (blue 
dot) compared to that in the SC state. Furthermore 
the resonance shifts to a lower energy (uj = I.6A0) and 
its width becomes broader, together with an intensity 
gain below u 1.4Ao- The behavior of x (Q: w ) is 
thus qualitatively consistent with the experimental ob- 
servation (see Fig. 1 in Ref. HH). Figure [TJb) plots 
X (q, uJres) in the SC state. It strongly peaks at q = Q 
where the spectrum is two-fold symmetric and shows 
a slight transverse-elongation, similar to that measured 
in Ref. [42]. In the mixed state, the spin excitation re- 
mains commensurate and x (q;W res ) does not change 
much except for that at q = Q, where the intensity 
is strongly reduced [Fig. HJc)]. Figure [TJd) shows the 
field-induced intensity change defined as A% (q, ui res ) = 
X (q,w res )| mixed— state X (q, w res )\sc- state- We can 
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FIG. 1: (color online) The s± case, (a) Normalized \ (Q,cj) 
as a function of the reduced energy (b) and (c) show 

X (q, w r es) in the SC and mixed states, respectively, (d) The 
field- induced intensity change (q,u; res ). (e) and (f) show 

X (q, lj) along the (q^O) direction at T — 0, for cj = 0.6Ao 
and lj = Lo rea , respectively. In (a), (e) and (f), I (blue dot), II 
(black solid) and III (red dash dot) represent the mixed, SC 
and normal states, respectively. 



see that the intensity is lowered around q = Q and 
the spectrum also elongates transversely. At T = and 
lo = O.6A0 (near the spin-gap energy), constant-w scan 
along the (qz,0) direction [Fig. [TJe)] suggests that the 
scattering in the SC state is very weak, with a small peak 
around q = Q, consistent with the presence of a spin gap. 
The magnetic field induces scattering in the spin gap and 
enhances the small peak around q = Q since there are 
in-gap Andreev bound states emerging close to the vor- 
tex centers^ In contrast, the imposition of the magnetic 
field at T = partially suppresses the resonance intensity 
at oj = UJ res [Fig. [T{f)]. 

In the s-wave case, we set U and Jh to be 4.1 and 0.5U, 
respectively, for the same reason as mentioned above. If 
we still choose U = 2 and Jh = 0.2U, no resonance will 
appear below u rj 5A in the SC state. On the other 
hand, if we choose U = 4.1 and Jjj = 0.5U for the s± 
pairing, magnetic order will appear in the SC state and 
coexist with superconductivity. However for the same U 
and Jh in the s-wave case, no magnetic order exists in 



FIG. 2: (color online) The s-wave case, (a) to (f) are the same 
as Figs. QJa) to[]Jf), respectively. The inset in (d) shows the 
spectrum close to q = Q. 




FIG. 3: (color online) Af(i—j, oj res ) for the s± (a) and s-wave 
(b) cases. 



the SC state since it cannot coexist with s-wave supercon- 
ductivity^ and this has been verified by our numerical 
calculation. Therefore we choose different U and Jh for 
the two pairing symmetries. The behavior of \ (Q,w) 
[Fig. [H[a)] shares some similarities with that in the s± 
case. The excitation is gapless in the normal state while 
a spin gap below lo rj 1.4Ao and a resonance at 10 = to res 
show up in the SC state. Upon applying the magnetic 
field, the intensity of the resonance is suppressed and its 
energy is lowered to lo = 1.6Aq. Similarly the resonance 
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becomes broad and below u> — I.6A0, the intensity in- 
creases due to the emergence of in-gap Andreev bound 
states near the vortex cores.— \ (q, w res ) in the SC state 
[Fig. [2fb)] also peaks at q = Q, with four- fold symme- 
try in contrast to the two-fold one in the s± case. In 
the mixed state, the peak in \ (ij^Ves) at q = Q is 
suppressed and the excitation area enlarges slightly [Fig. 
EJc)]. The field- induced intensity change in Fig. EJd) 
shows that, close to q = Q, the intensity is suppressed 
by the application of the magnetic field. However, further 
away from q = Q, there are intensity gains in the red and 
yellow areas. This behavior is different from that in the 
s± case where there are no intensity-gain areas around 
q = Q as shown in Fig. [Tfd) . Thus the distinct char- 
acteristics in the field-induced intensity change can be 
used to distinguish the s± and s-wave symmetries. In 
the spin gap, the constant-w scan along the (q x , 0) direc- 
tion at T = and uj = O.6A0 [Fig. E£e)] is qualitatively 
the same as that for the s± pairing, where the magnetic 
field enhances the small peak around q = Q due to the 
Andreev bound states. Contrastingly, at the 
magnetic field suppresses the resonance at q = Q but 
enhances the intensity from c^/n = 0.7 to 0.9. 

The behaviors of x (Q? w ) shown in Figs. [IJa) and 
[Ha) can be explained as follows: In the SC state, the 
resonance energy is approximately at u> res = |A(k)| + 
A(k + Q)|, where k and k + Q are the momenta on the 
hole- and electron-pockets, respectively Since the aver- 
age SC OP in the mixed state is reduced to 0.95Ao, thus 
the resonance shifts to a lower energy. Furthermore the 
application of the magnetic field leads to the emergence 
of in-gap Andreev bound states which induces scattering 
in the spin gap, thus the spectral weight is transferred 
from the resonance to that in the spin gap, leading to a 
weaker and broadened resonance. 

On the other hand, to understand the different features 
of Ax (q, UJres) m Figs. (Hd) and EJd), we study the 
properties of Af(i — j, u) res ) defined as A/(i — j, u) res ) = 

f\i> ~ 3i Wresjlmixed- state — f(i — j,^res)\sC- state- For 

the s± pairing, the mixed-state effect influences only the 
short-range spin-spin correlations [Fig. EJa)] and the 
spectrum in Fig. HJd) can be reproduced by considering 
only Af(i - j, u res ) in the range \i x - j x \, \i y - j y \ < 4. 
For the s-wave case, in contrast, the mixed-state effect 



extends to much larger spatial areas and all the points 
in Fig. [3jb) have to be included when performing the 
Fourier transform in order to match the spectrum in Fig. 
[5Jd). Thus if the real-space spin-spin correlation func- 
tion can be measured experimentally, Fig. [3] can be used 
as another method to differentiate the s± and s-wave 
pairings. 

Summary. — In summary, we have systematically stud- 
ied the mixed-state effect on the low-energy spin dynam- 
ics in optimally-doped iron pnictide superconductors, for 
the s± and s-wave symmetries, respectively. The spin ex- 
citation spectra at q = Q are similar in these two cases, 
no matter it is in the normal, SC or mixed state and the 
evolution of the resonance with the magnetic field are 
both consistent with the experimental observation, that 
is, at zero field, a spin gap and a resonance show up in the 
SC state while in the mixed state, due to the existence 
of Andreev bound states and the suppressed SC OP, the 
spin gap is filled up and the resonance shifts to a lower 
energy and becomes broad, together with a reduced in- 
tensity However, the field-induced intensity change at 
the resonance energy shows different behaviors between 
the two symmetries. For the s± pairing, the intensity 
of the spin excitation is suppressed around q = Q in 
the mixed state while for the the s— wave case, there are 
intensity- gain areas in the vicinity of q = Q. Further- 
more, the mixed-state effect affects only the short-range 
spin-spin correlations in the s± case while it extends to 
much larger spatial areas in the s— wave one. We thus 
propose to measure the field-induced intensity change in 
momentum space and if possible, in real space, to distin- 
guish the s± and s-wave symmetries. 
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